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Inclusion of the Third Zonal Harmonic in an 
Accurate Reference Orbit of an Artificial Satellite 

John P. Vinti 1 

(November 17, 1965) 

The previous paper gave the general theory and physical principles involved in so modifying the 
author's spheroidal potential for an oblate planet as to permit exact inclusion of the effects of the third 
zonal harmonic of the planet's gravitational field. The present paper carries out the computational 
details necessary to derive the resulting orbit, which now corresponds to a potential fitted exactly 
through the third zonal harmonic and to about two-thirds of the fourth. 

The accuracy of the orbit itself, as a solution for the given potential, depends on the accuracy of 
solution of a certain cubic equation. The paper works out this solution through terms of the third 
order in A, the coefficient of the second zonal harmonic 1 , but its accuracy, and thus that of the secular 
terms, may be increased at will. Periodic terms arc carried through the second order, but their ac- 
curacy may also be increased. 

An obvious advantage of accounting for ,/ :! in this way is the absence of small denominators in 
e or sin I that occur in a perturbation theory. Another is the resulting increase in accuracy, through 
terms in ,/':■, of the long-periodic third harmonic terms. 

Key Words: Third zonal harmonic, increased accuracy, spheroidal reference orbit, artificial 
satellite. 



1 . Introduction 

The previous paper 2 [V1966| developed the theory and physical principles for inclusion oi 
the third zonal harmonic of a planet's gravitational potential in an accurate reference orbit of an 
artificial satellite. The aim of the present paper is to work out the resulting generalized reference 
orbit in such detail as to permit practical calculation. 

To begin, Jet p, 77, if be oblate spheroidal coordinates, satisfying the equations 

x + iy=r cos 6 exp ^ = [(p 2 + c 2 )( 1 — t? 2 )] ,/2 exp up (1.1) 

z=r sin 0=— S-r-prj. (1.2) 

Here x, y, z are rectangular coordinates of a satellite in a Cartesian frame, with origin at the 
center of mass of an oblate planet, Oz pointing along the polar axis, and Ox pointing toward its 
vernal equinox. Also r, 0, <p are the planetocentric distance, latitude, and right ascension. 
If r e is the equatorial radius and if 

c 2 =Ai/ 2 (l-^ 2 ./") (2) 



8 = -h'J 2 \J:u (3) 



then the gravitational potential 



K = -^(p 2 + cV)-'(p + T7S) (4) 

leads to separability of the problem of satellite motion [VI966J. Here 8 — 7 km for the earth and 



"Research supported by the National Aeronautical and Space Administration, Washington, D.C. 

1 Present Adress: Department of Mathematics, North Carolina State University, Raleigh, North Carolina. 

' Since there will he so many references to my own papers. I ahhreviate my name to the initial V in the body of the paper. 

17 

796-485 O-66— 2 



the potential (4) leads to a fit of 



l-JP JnP„(sine) 



(5) 



exactly through the third zonal harmonic. It also gives about two-thirds of the fourth harmonic 
and negligible values, of order J\ or smaller, for the higher harmonics. 

By V1959, if a.\ is the energy, a 3 the z-component of angular momentum, and a 2 the separation 
constant, the generalized momenta p p , p^, and p$ are given by 



p* = a 3 (6.1), p p = ±(p 2 + c 2 )- 1 F(p) 1 / 2 (6.2), p 1 = ±a-7 ] '^G(r } ) 

Here F(p) and G(tj) are the quartics 

F{p) = c 2 al + (p 2 + <?){- «1 + 2p,p + 2a lP 2 ) 
G(ij) = - a\ + (1 - r] 2 )(al + 2p^8 + 2a 1 c 2 T/ 2 ). 

The Hamilton-Jacobi function JF(p, 17, </?) is then 



1/2 



(6.3) 

(7.1) 

(7.2) 



W= [p*dy+ [p P dp+ J pvdr ] =az<p + j P ±{p 2 -^ c i )- 1 F(py' 2 dp + [^ ±{l-rj 1 Y'G{r ] )^dr). (8) 

The quartic F{p) has two physical zeros, p\ and p 2 , and the quartic G{y]) has two physical zeros, 
j)\ and 170- The possible motions are thus confined to the toroidal region 



0<p!^p^p 2 
— 1 ^ 771 ^ J] ^TJo^ 1. 
If /3i, /3 2 , an d /^ are constants, the orbit is then given by 



5a 1 da>> dati 



With pi and 171 as lower limits of the integrals in (8), these become 

t + p^Ri + c'N! 

if = ft 3— c 2 a : ^3+ a-.iNr 

Here the p-integrals are given by 

Rr= f P ±p 2 F-" 2 dp, R 2 = [ P ±F-^dp, R 3 = f^^ + c 2 )- 1 / 1 " 1 / 2 
Jpi Jpi Jp\ 

and the ^-integrals by 



(9.1) 

(9.2) 

(10) 



dp 



A'; 



J'h 



(11.1) 

(11.2) 
(11.3) 

(12) 

(12.1) 

(12.2) 
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/v, 



= \±(l-rf)- l Gr 1 l*dri' 



(12.3) 



With T7 1 as lower limit in the 77-integrals, rather than as before, we shall find that in the 
Keplerian limit U-i^ 0, J 3 — > 0) the constant — fi\ will again reduce to the time of passage through 
perigee. The constant fi' 2 , however, will not reduce to a>, the argument of perigee, but to co + 7r/2. 
Similarly /% will not reduce to ft, the right ascension of the ascending node, but to ft — $77 sgn a 3 . 

The p-integrals have already been worked out [V1961b|. The new ^-integrals, however, 
will contain 8, involving / 3 , and will thus have to be reevaluated. The task of the present paper 
is thus to reevaluate these 77-integrals TVi, /V 2 , and N3, insert them into (11), invert (11.1) and (11.2) 
to find p and r\ as functions of time, and then insert the results into (11.3) to find if as a function of 
time. We shall then have the new reference orbit, corresponding to the generalized potential 
(4) that allows for exact inclusion of the third harmonic. 

2. Factoring the Quartics 

In order to evaluate the integrals, we first have to factor the quartics F(p) and (7(7?). The 
factoring of F(p) was essentially done in V 1961b and the only new task that we have in connection 
with it is to relate associated quantities to the new orbital elements, which will depend on the 
new factoring of 6(77). To factor (7(7?) note that by (7.2) 



(7(7/) = (XTy — OL^+ 2fJLT\?)+ {2tX\C~ — ocr^j] 1 — 2/x7/ :} 8 — 2aic 2 7) 4 



(13) 



By (6.3) and (13), if tin 



rbit 



reaches r\ = 0, we should have of — o§i^0. 11 



d tl 



i\d 



harmonics are all put equal to zero, we should have 8=2 r ( J 2 ] | ,/.$ | =0 and by (1.2) a motion z=0 
in the equatorial plane would be possible, with 17 remaining zero and a| — al ^ 0. 

As soon as we insert the third harmonic, however, there arises a doubt that 77 = will always 
be reached. Accordingly, we must leave the sign of ofi — a$ open. For bounded orbits we now 
have ol\ < and the table 



-1 


+ 1 

+ 00 



Qv) 



a l ~~ a i- sig" either 
- a'i ^ 



There are thus two non-physical zeros, between — oc and — 1 and between + 1 and +°°, and two 
physical zeros, between —1 and -hi. Call these latter 771 and t/ () , with 771 = t/ . If 8 = 0, G{tj) 
becomes a quadratic in rf and 7/1 = — 7/0, where 170= 0, the equality sign holding when d\ — a§ = 0. 
To see what happens when 8 7^ 0, make a plot of G(r)) versus 77, noting that G"(0) = 2/x8 > 0. 
In order to avoid more than two zeros between 07=—! and 17 = +!, we must then have 
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if that and 

a\ — «1>0 i7i<0<77o — T7oi?i>0 

a 2 2 — al = 0=T7i<i7o — t7oT7i = 

al~a 2 3 <0 0<7)i^T7o -t?oT7i<0 

The paper V1961b indicated the following form for factoring G(r}): 

G( V ) = (a\ - a%) (l - ^j (l - ij (1 + Cut - C 2 r? 2 ). ( 14) 

Here, since 171 ^ 17 ^ 170, we have 

(17-171) (t?o-t?)^0. (15) 

For af — a| > 0, we have —r) {) 7] l > and division gives 

For ai — a$ < 0, we have 77017 1 > and division gives 

K)K)<°- 

In either case, since £(t)) is 0, we find that the factor 1 + C\Tj — C217 2 > 0, where C\ and C2 will 
both be small compared to unity. 

On comparing coefficients in (13) and (14), with 

(18) 
(19) 
(20) 



(21.1) 

(21.2) 



2P+C l S = — S (22.1) 

go 

l + C i S-2C l P = Pu + c2lao S (22.2) 

go 

C I + 2C 1 P = -S (22.3) 

9o 
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go = fJL'Hal— a%) 








P = 2 (l7o-r-T7i) 








5 = — 170171, 






and with the relations from V 1961b, 


viz, 






i V* 

(Iq =— 2 , 

«1 


_/ 2a 1 oty/ 2 
e = IH —J , 

\ ^ / 

a| = MPo, 


Po = a (H 


-ej) 



we nave 



c 2 

c 2 =—s- 

cioqo 



(22.4) 



Along with three constants, to be derived later from j8i, /3!>, and /?.,, we shall choose as orbital 
elements 

a= ±(pi+p2) (23.1) 

e = ^ (23.2) 

p = a(l-e 2 ) (23.3) 

as in V1961b, and 

S=-i7oi?i. (23.4) 

[f8 = 0, T}i becomes — 170, and S reduces to the old rjf ) = sin 2 / of VI 961 b. Note that now, however, 
the new orbital element S can sometimes become negative, viz, for orbits sufficiently close to 
equatorial. Note that for a polar orbit 770 = -\- 1 and rj\ = — 1, so that S then equals unity. In any 
case, since — 1 ^ 171 ^ 170 = 1, we have 

S 2 =il. (24) 

Our next task is to express the factorization of the quartics F(p) and G(r)) completely in terms 
of the elements a, e, and the new S. To do so, we must first express Ci, C2, and P in terms of 
c, 8, ao, po, and S; that is, we must eliminate q^. 

On inserting (22.4) into (22.3) and solving the resulting equation simultaneously with (22.1), 
we find 

c ,.< 1 -*- S .)"»sh-*-s) ,25.1, 

P=(l — — S'V'-SU-S). (25.2) 

\ aoqo I qo 

If we now insert (25.1), (25.2), and (22.4) into (22.2) and call 

**S-u, (26) 

we find 



wfa-nfl-JL) 



I =1+ _cl (1 _ S)+ W^ ,>" »W . 27) 

// a />o /I r" nV 



For 8 = 0, S becomes sin 2 /, w becomes sin 2 //sin 2 i , where 



io - cos" 1 — , (28) 
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from V1961b, and (27) reduces to 



c 2 



sin 2 i = sin 2 / H sin 2 / cos 2 /, (29) 

CIqPo 

in agreement with (4.8) of V1961b. Equation (27) is a cubic in u~\ but we may solve it by succes- 
sive approximations. To do so, note first that 

-=1 + -^(1-S) + 6>(S 2 ). (30) 

u aopo 

If we insert (30) into the right side of (27), the error will be of order 8 4 , i.e., of order ]\. Thus 

-=1+ — (l-S) + €, (31.1) 

where 

'**)\l-S)(l--±-s) 

^ 2 V %° ' + OUi). (31.2) 

1 + — (1-2S) 

«oPo 

(Insertion of (31) in the right side of (27) would reduce the error to 0(/£)> but such accuracy is not 
needed.) Equation (31) frees u from q and thus also Ci, C 2 , and P, since by (22.4), (25), and (26) 

c 2 
C 2 = — u (32.1) 

a po ' 

C, = (l — Su) — ttfl —u) (32.2) 

P=(l- — SuY'-uil-S). (32.3) 

V a po / Po 

These quantities Ci, C 2 , and P are then correct through order J|. 
Since |S| = 1, we find from (31) that e i^ and ^ _1 i^ 1. Thus 

0<w^l (33.1) 

and then, by (32) 

C 2 >0, Ci>0, Pi=0. (33.2) 

In order to obtain all the required quantities as functions of a, e, and S, we must now summarize 
the factoring of P(p), given in V1961b. We have 

F( P ) = -2a,(p- Pl )(p 2 -p)(p 2 + Ap + B) (34.1) 

« = i(pi+p 2 ), e-^-^ 1 , p-a(l-e 2 ). (34.2) 

p 2 + pi ; 
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Then, with q defined by (18) and a and p by (21), comparison of coefficients in (34.1) and (7.1) 
gives 



2a-A = 2a 



B + ap — 2 A a = c 2 + aopo 



2aB-Aap = 2a c 2 



Bap = (i()P()( :1 (i()lpo- 



(35.1) 
(35.2) 
(35.3) 
(35.4) 



We now regard a, p, S, c, and 8 as known and have to solve the eight equations in (35), (32), and 
(26) for the eight unknowns a () , po, ao, A, B, Ci, C 2 , and P. Here a is an auxiliary parameter given 
by (31). 

From (35.4), (35.2), (26), and (31), we find 



Hap = c 2 S 



4^5 2 (l-,S)(l- 



B + ap — 2Aa — c 2 S -f 



Po 



(loPo 



+— (1- 

dnpn 



■2S)J 



OUt) 



Now since /4 and i? are both of order ,/ 2 , we have from (35.1) and (35.2) 

c 2 c 2 

=—+(KJl) 

aopo ap 



«o = 2( 1+ .ci+M_d_JL') f0(J|) 

Po p V ap p a ap 



(36) 



(37.1) 



(37.2) 



Because of the factor 8 2 , of order J 2 , insertion of (37) into (36) will not change the order of the error. 
Thus 



#up = c 2 S \b + a/> - 2^a - c 2 S + — (1 - S)8 2 { 1 + — (3S - 1 ) + A (- 



ap 



) a J (ip)_ 



+ 0(71). (38) 



This is an equation i 
of (35.1) and (35.3). 


nvolving only A 
It is 


, B, and "known" quantities. 


Another 


follows 


by 


comb 


ination 








2aB — Aap = 


(2a-^)c 2 . 












(39) 


The 


solution of (38) 


and (39) for A and B is 


















A 


— 2ac 2 (ap — c 


•S)(l s) + 8a2< 
p 




(3S- 


-2)JS(1 


-S) 






(40.1) 


(ap — c 2 ){ap- 


c 2 S) + 4a 2 c 2 S4 


4c 2 

— ;-8 2 (3ap- 
P 


-4a 2 


-c 2 )S(l 


-S) 








B = (2a)- 1 (op- 


-c 2 )/4 + c 2 . 












(40.2) 








23 

















For a polar orbit, for which S=l, we have from (40) and (35) 

A = B = c 2 a = a Po = p. (41) 

For such a polar orbit, if 6 = 0, we then have 

p = a = ao = p = po, e = 0. (42) 

Equations (40) give us A and B in terms of known quantities. To find the others, note that 
do then follows from (35.1), p from (35.2), q from (26) and (31), and C2, Ci, and P from (32). With 
use of these equations, the errors arising from the approximate solution of the cubic (27) are of 
order J\ for all the unknowns except q . For q the error is of order J|. The "known" quantities 

are c, 8, a, />, and S. Here c = r e J} 2 12 f 1 — -jJiJ^ 3 ), $ = 2 r e ./ 2 " 1 | t / 3 |, a is a constant "semi-major 

axis," p a constant "semi-latus rectum," and S a constant quantity analogous to the squared sine 
of the inclination for inclinations not too close to zero. For inclinations close to zero, however, 
S may be negative for 5 7^ 0. 

3. The 77-lntegrals 

From (14), (18), (19), and (20), we find 

G( V )= f xqoS-HS + 2Pr ] -r ] 2 )(l + C i r } -C27 1 2 ). (43) 

Now 

S + 2Pr/-ry 2 ^ J P 2 + 5-(7 ? - J P) 2 . (44) 

The physical zeros 771 and 770 of G(r}) are the zeros of (44), viz, 

Vo = P + iP 2+Syl 2 (45.1) 

Vl = P-(P* + S)^. (45.2) 

For 8 7^ 0, S may be negative, with smallest possible value S = — P 2 , corresponding to 

from VI 966. 
With 



Q = (P* + S)" s &0, 



(47) 



we may now introduce a variable t//, such that 

7j = P+ Q sin ifj. (48) 

If 1// continually increases, so that 1// > always, 17 will then vary between the limits rji and 170 and 
the double sign in eqs (12) for the ^-integrals will disappear, with 1// as integration variable. Thus 
\\) is a uniformizing variable. 
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Insert ion of (48) into (43) then gives 

G( v ) = iJLqoQ 2 S Hl + C ]V -C 2V 2 ) cos 2 i//. (49) 

and the factor ±G~ l/2 dr] in the integrals (12) becomes 

± G- X i 2 dy] = (-^-gj * (1 + Ciri - C> v 2 ) l ' 2 d0 (50. 1 ) 

= oq x u x H\ + Cir\ - C 2 t? 2 )- 1/2 # (50.2) 

by (26) and (21.2). 

4. The Integral /V2 

By (12.2), (50), and (48) we then find 



7_ 

"2 



a 2 u- 1 l 2 N 2 = I jl + Cr)-^ 2 )- 1 ^. 



Since \ = 0(J¥ 2 ), it follows that |X± V\ 2 - 1 1 = 0( 1) and 



(51) 



Here C\ and C 2 are both positive, of order J-,. To evaluate (51), we plaee 

I, = ,;CV 2 r, = (>(./,), X - I i&Ci* = ()(.}.)% (52) 

so that 

1 + Ci-q - C 2 r} 2 = 1 - 2XA + ft 2 . (53) 



| A | ^smaller of |X±Vx 2 -l~|. (54) 

Thus a Legendre expansion is legitimate, so that 

(l + C 1 r ? -C 2 r/ 2 )- 1 / 2 = (l-2X/r-h/ ? 2 )- 1 / 2 - J h n P n (k)= J i»C$ 2 r)»Pn(k). (55) 

N = /) = () 

Then (51) and (55) lead to 

a 2 u-^N 2 = J F„(Ci,C 2 ) [V#, (56J) 

where 

F»(Ci,C«) = i"Qi 2 P„(h tC^" 1 ' 2 ). (56.2) 



are real polynomials in Ci and C 2 . Through F 8 they are 

F.= l, F,— |Ci, F 2 =§C?+fc 2 , F 3 =-^C?-|CxC 2 , 
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at — ^ rn . 315 , 4 ^ 105 ^p ± An F —— -=- r? — ^5 r^C — ^^- T3 r2 — — r T3 

n 1024 Cl ^256 Cl ° 2 ^ 64 ClC2 + 16 C2 ' ^ 7 2048 Cl 512 Cl ° 2 128 ClL2 32 Cl ° 2 ' 

_6435 _ 3003 ^ 6 r 3465 ^, 4 ^ 2 315 ^, 2 ^ 3 35_ 
^s- 12g c i + 2048 ClC2+ 1024 ClC2+ 128 Cl 2+ 128 Cy (57) 

Higher F's are of order J I or smaller. 
We now have to evaluate 

J ^ n dijj = J ^(P-Q S mx) n dx== j^ (P+0cosy)»rfy=/ n (u/), (58.1) 

~ 2 ~ 2 

where 

tt? = tt/2 + i//. (58.2) 

As in V1961b, p. 178, we now readily show that 

f n (w) = 7T- l wf n {7T) + g n (w), (59) 

where gVi(w) is an odd function of w, of period 2tt. From that same paper 



f n (ir) - f * (P ~ Q cos y)"dy= ttP"(1 -ef) M/2 Pn[(l - ei)" 1 ' 2 ], 
Jo 



(60.1) 



where P^M is a Legendre polynomial in x and where 

e x s £/P (60.2) 

< ei < 1 for S < (60.3) 

d > 1 for S > 0. (60.4) 

For S < 0, the factors (l-e?)" /2 and P n [(l - ei)- 1 ' 2 ] are both real, for all n. For S > 0, they are 
both real if n is even and both imaginary if n is odd. In any case their product is real. 
If we now introduce, as in Y1961b, p. 178, 

Rn(x) = x n P n (llx\ (61) 

a polynomial of degree [n/2] in x 2 , we find 

7T-%(n)=G n (P v Q) (62.1) 



G n (P, Q) ss P"R n (VT^\ (62.2) 

polynomials of degree n in P and (). Through rc = 8 they are given by 

Go- 1, d = P, G 2 = P 2 +|<?, Gs^P^ + Ip^? 2 , 
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G 4 = P* + 3P«<?H | Q\ G s = P->+ 5P 3 <?* + ^ PQ*, 

16 C« = 16P" + 120 P*(>- + 90P 2 *? 4 + 5(? K , 166' 7 = 16P 7 + 1 68PC? + 210P^ + tfPQe 

128 C 8 =128P 8 + 1972P«(? 2 + 3360P 4 <> 4 + l^OP-'C^ + SS^". (63) 

On inserting (58), (59), and (62) into (56.1), we find 

a 2W -" 2 A r 2 = (| + «|/W+ J P„(C. C 2 ) J c Hj • sin>, (64.1) 

where 

fi 2 = JF„(C, C t )G„(P,Q). (64.2) 

Remembering that G, C 2 , and / J are all of order / 2 and that Q is of order /§, we have from (64.2), 
(57), and (63), through order Ji: 



P,= 1 -|c,P + (| Cf +|c 2 ) (P 2 + |<7) +^ ci<? -| C,C2P^ 

-Iw+lw+lw+^CK!' 



Q 

Note that Z?> contains — 2C1/H 77 C\Q Z , proportional to J\IJ\ = O(M) , so that inclusion of ./. } in 

Lo 

the separable problem apparently leads to new second-order secular terms. We shall say more 

about this point later. 

Let us check Bo for 8 = 0, in which case it should reduce to 

fi2 = ^/^) = l + I^ + |^ + |^ + (^)V. • -, (66) 

given on p. 184 of V1961b. For 8 = 0, we have C\ = 0, / >= =0, r)\ = — 770, S = j]l = sin 2 /, and Q = y] {) 
= sin /. For Cj we then find from (22.4) 

^2 = — — m = ~> — 2 ^ < 67 ) 

a oPo g "<>/>(> cr> — «"} 



nth use of (18) and (21.2). But in the notation of V1961b, {di - (4)1(4 = sin 2 to, so that we find 

Cz = _£!__j|_ (68.1) 

tfoPo sin - io 

= iff«, (68.2) 
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from p. 173 of V1961b. On inserting these values into (65) and using q = 170/172, we find 

This agrees with (66), so that B 2 reduces correctly for 6 = 0. 

We now have to find the periodic terms in (64.1), i.e., the coefficients c n j- Since F 5 (ci, c 2 ) 

is of order 71, we need take n only through the value 4. By expressing (PH- Q cos y) A , A; = 0, 1, 2, 3, 4, 

f w 
as trigonometric polynomials, evaluating (P — Q cos yfdy, picking out the periodic terms, and 

Jo 

keeping results only through order 71, we find from (56.1) and (57) 

a 2 u- 1 ' 2 N 2 = (|+i//W + (P-^Cx(Ac*Q sin w 

+ i [(4C 2 + 3C1)<? + 3CK) 4 ] sin 2w + ^ Q> sin 3w 

oZ 16 

3 

+ -37 CiQ 4 sin 4m; -f periodic terms of order /§. (70) 
zoo 



In terms of 1// = m; — 7t/2, this becomes 



a 2 tt- 1/2 #2=U + <M B 2 



(=+♦)* 



p + j|c,0,)+|c,<? 



COS \\) 



- ^ [(4C, + 3Q0 2 + 3C1CH I sin 2«fr - ^ O 3 cos 3 >// 

+ ^7 CiQ* sin 4»/; + periodic terms of 0(71). (71) 
p or § = 0, we find that the periodic terms in (71) agree with those on p. 184 of V1961b, as they should. 



5. The Integral N t 

By (12.1), (50), and (48), we have 



a 2 M-" 2 /V, = f * ^ (1 + Ct; - C 2 t) 2 )- 1/2 tj 2 #. 



(72) 



Then, from (55), (56.2), and (58.1) 

a.u-^N, = J F„(Ci, C 2 )/„ +2 (w). (73) 

n=0 

As before, 

/ B+2 (u>) = wG n+2 (P, Q) + gn +i (w), (74) 
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where g,i+2(w) is odd in W, with period 2tt. Thus 

aar^Ni = (f + $)B\ + J FJL& 9 C 2 )#„ +2 (k;), (75.1) 

where 

flj - J F„(d, C 2 ) G„ +2 (/', £>)• (75.2) 

Since the integral TV] has a factor c 2 = 0(j2) in the kinetic equations (11), we need each F n G n +2 
correct only through order J|, in order to obtain secular terms correct through order J\. From 

(75.2), (57), and (63), we then find 

fii= i@2 + p2_| Cl P22 + | C2 p202 + A (4C 2 + 3C?)C? 

4 z o4 

-i ClC ^ + 2f6 (6 ^ + IWW+^CW + OW). (76) 



For 6 = 0, this reduces to 



3 15 175 



or 



W-*+£^l|* + 55S* + . • • (77) 

in agreement with the value of B\ on p. 184 of V1961b. 

To find the periodic terms in (75.1), we proceed just as for /V2, obtaining 

a 2 «-" 2 /V, = (| + <fr) B[ + (ziV + l w) cos 4> 

- 1 (2<3 2 + CiQ*) sin 2x\, - ^7 C1O 8 cos 3i|/ 
o 24 

+ Z7 CzC? 4 sin 4i// + periodic terms of ()(fi). (78) 
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For 6 = 0, the periodic terms agree with those appearing in (6.38) of V1961b. 

6. The Integral /V» 

By (12.3), (50), and (48), we have 

a2// -l/2/V, = ( l_^ ) -. (l4 _ Ci7? _ C2T? 2 ) -l/2^ (?9) 

~ 2 

Then, from (55) and (56.6), 

a 2 u- 1 l 2 N 3 =^F n (CuC 2 )L n (80.1) 

n = 

29 



7)" 1 

7) I-17 



n" (—1)" "~ l 
1 + 7J 1 + 7, 2,/ "> ^ * 



Addition of (83) and (84) then gives 



reducing to the two cases n =odd and /? = even, viz, 



a 2<T 



~2(r+l ^ 

1 - rj 1 = 1 - T] 2 



The lemmas (82) then follow directly from (80.2) and (86). 

With application of (58.1), (59), and (62), we then find from them 



a— 1 a- 



(80.2) 



(81.2) 



Ln= /"* (l-77 2 )77«#. 

In particular 

/>o = J (1 — T7-)-^/i// (81.1) U = J (l-T) 2 )- 1 ^^. 
To evaluate the L,'s we first prove the lemmas 

L 2 . + i = Li-2 P *? 2r+1 # (cr^l). (82.2) 

To prove them, note first that 



!_„ ,— "S^ (83) 



(84) 



» ^H J (—1)" »- ] 

T J3 - + T^ = 7-^ + i rr^-S [l + (-l)" + '"1 17*, (85) 

1 — 7) 1+17 1 — 77 1 + 7) jgj 



7)^" « — J 



T) 2<r+1 n 

- £ t,^, (o- s l). , 86 .2) 



L 2<r = U- w £ 6', 7 (P. C) - X ** M ' °" - l (87-D 

r=0 r=0 

Z.2, r+ i = L 1 - w < 5'G 2 r + i(P, (fi-^gn-iiw), °" = * > (87.2) 

T=0 T=0 

where w = ijf + 7r/2. Insertion of (87) into (80.1) then gives 
a 2 u- l l 2 N 3 = L J Faor + Lj J F> 1(r+l + B*w 

(7 = (7 = 

- J f 2,r J ^W) - J ^*H"1 2 ^r + l(w), (88.1) 

(7=1 7=0 (7=1 7=0 
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B 3 - - J F 2(r 2 C* - £ F 2 „ +1 X ftr+i. (88.2) 

rr = l r=0 a = l r=0 

To simplify the first terms in (88.1), note that by (56.2) and (52), if 

h = iO/ 2 , (89) 

the value of h for r] = 1, we obtain 

F 2(T = h^P^W (90.1) F 2(r+1 = Ar +1 i>2<r+i(X). ( 90.2) 

To evaluate the sums VF 2(r and VF^+i, note that 
o o 

a-2\ho+hi)-n*= J h*rp 2 ^k)+ J / l( r +i ^,r + ,(X) ( 9i.i) 



and, since P 2 <r(— A.) = P 2o -(A.) and F 2(r +i(— X) = — F 2(r+ ,(A), that 

(1 + 2U +A§)- 1 / 2 = J ft8"fto(X)- J fcj* +1 JWX). (91.2) 



Solving eqs (91) for each sum, with use of (90), we find 

2 J F 2(r = (l-2X/io-r-/i§)- 1/2 + (l-f2XAo+A8)- 1/2 (92.1) 

o 

2 J]F 2(r+1 =(l-2X//o + /^ ) )- ,/2 -(l + 2X/ ? o-f/4)- 1/2 . (92.2) 

o 

Since /i = /?<> when 17 = 1, we have from (53) 

1 T 2kh -f /if) = 1 ± C, - C 2 . (93) 

Then from (92) and (93) 

U J) F 2(r + L, 2 F 2(r+1 = i L [(l + Ci-Cj-w + a-Ci-C*)- 1 !*} 

z 

-f|l,[(l +C 1 -C 2 )-^-(l-C 1 -C 2 )- 1 / 2 ]. (94) 



For Z? 3 , the coefficient of a secular term, we find from (88.2), (57), and (63), carrying terms 
through order /|, 
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cs(i+|<?'+l<? l +^0')-(^ci+|cift)i' 



■jffcqfr+f 



+ £ C, Q P + 1 P<? + 0(JI). (95) 



Here the term — — C\ is proportional to J\l J\ = 0{J\) , so that we apparently again find a second- 



3 

8 
order secular term produced by inclusion of J3 in the separable problem. 

To check (95) for 8 = 0, put Ci = P = 0, C2^T7~ 2 , and Q = r)o. It becomes 



^3=|^+|^(l + |^) + ^^(l + |^+|Vo) 

+ ^ r, 2 -« (1 + 2 ^0 + g rt + ^ tj§) + 0(71), (95. 1) 

agreeing with (6.65) of V1961b through terms of order J\. 

To evaluate the periodic terms in (88.1) through order J|, note that all coefficients F n are of 
order J\ or higher for n ^ 5. Thus we need only go, gi, and #2, where g n is the periodic part of 

fn= j™ (P-Q cos y)*dy. 



Then 

Q 2 

go = 0, gi = Q sin w, g 2 = — 2P() sin w + — sin 2w. (%) 

From (88.1) the relevant periodic terms are then — F 3 gi and —Fig*. With the aid of (96) and (57), 
their sum then becomes, through order Jf, 

3 3 

Periodic terms = — - C\C 2 Q sin w— — C\Q 2 sin 2w (97.1) 

= -|c,C 2 <?cos * + ^ Cf<? 2 sin 2i//. (97.2) 

We come next to the evaluation of L<, and L\. From (81), we have by expansion in partial 
fractions 

2L = M' 1 + M' i 

2L 1 = M[-M' 2 , (98.1) 



whe 



r;=f* (i-tj)- 1 ^ (99.i) M^=f* (i+7))-'#. 

J -ir/2 J -tt/2 



With the aid of (48) and the definitions 



e ^Y3p (loo.i) e s s y+p' (100 - 2 > 
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we find 

A/; = (l-/V (l-e 2 sin x)' l dx (101.1) 

J -tt/2 

M! 2 = (l + P) 1 {l + e*sinx)- l dx. (101.2) 

J 77/2 

Here e-z t== 1 and ez= 1, as we shall now prove. Since |S| ^ 1, we have by (31) that u ~ 1 and then 
by (32.3) that 

p=«i-S) (102.1) ^(i-^-s^A, (102>2) 

where 0< { < < 1. Now by (47) 

( 1 - Pf -Q 1 =\-2P + P 1 -P 1 -S=l- 2P-S = (1 -S)(l - 20, 

by (102). But 0<£<< 1, so that l-2£>0and 1 -S ^ 0, so that (1 - Pf - Q 1 ^ and 1 -P^Q, 
since neither \ — P nor (J is negative. Thus e 2 = l. Also (1 + P) 2 — Q 2 = 1 + IP — S ^ 0, since 
P^OandS^ 1. 
Thus l+P^Q and e 3 ^ 1. 

If in (101.1) we place x = y— 7r/2 and in (101.2) we place # = y-r-7r/2, we now find 

(103.1) 



M[ = (l-P) 1 I (l + e 2 cosy)- 1 rfy 

Jo 

/>-ff/2 

(l + e 3 cos y)-^y. (103.2) 



Af.: = d + P)- 



To evaluate the M"s we define a function E' e (y) by 

cos E>= ^ + C ° S y sin£,^ (1 r/ )1/2siny («<1). (104) 

1 + e cos y 1 + e cos y 

Then when y=/i7r, we find sin E' e — and cos £ , g = (—l) n . We may therefore impose the further 
requirement that E' e {y) — y whenever y=Ai7r. The relation between E' e and y is then the same as 
that between the eccentric and true anomalies, respectively, of an ellipse with eccentricity e. 
We next prove that 



J " (1 + e cos y)- l dy=(l-e 1 ) " 2 E' e (w). 



(105) 



To prove it, we first derive from (104) the relations 



(1-e 2 ) 1 / 2 sinE" 

sin y= — = — (106) 

1 — e cos A 



a 2 



— — 1 — ecosE'. <im\ 

1 -re cos y UU/j 

J3 
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Differentiation of (107) gives 



(1 — e 2 ) sin y dy 
(1 -f e cos y) 2 



sin E'dE'. 



Division of (108) by (107) then shows that 

sin y dy sin E'dE' 



1 + e cos y 1 — e cos 2? 



-=(l- e 2)-i/2 s ' my dE', 



by (106). Thus 



dy 



\-\-e cos y 



(1-e 2 )" 1 / 2 dE\ 



(108) 



(109) 



(110) 



from which theorem (105) follows at once, since E' = when y=0. Then, since ^ (0) =0 and 
E' e (—7r) =—tt, we find on applying (105) to (103), 



M\ = (1 - P)~ l (1 - ef) - 1/2 £ 2 (* + ir/2) 

WJ=(l + f)-l(l-«j)-^ [^(l//- 77/2) +77]. 

From (111), (100), and () 2 = P 2 + S, it then follows that 



(111.1) 
(111.2) 



M[ = (1 -S -2P)- l l 2 E! 2 (ifj + 77/2) 
M.;=(l-S + 2P)- 1 / 2 [^(i//-77/2) + 77]. 
By (112), (98), and (102), we then find 

77/2 



U-S) 1/2 Z,<r 

(l-S^L,: 



Vl + 2£ 

-77/2 

Vl + 2£ 



Xo 



+ Xi. 



(112.1) 
(112.2) 

(113.1) 
(113.2) 



9 _ ^(i// + 77/2) E^-ir/2) 



2Xi 



Vl-2£ Vl+2£ 

.^W + 7r/2) EM-tt/2) 



Vl-2£ Vl+2£ 

From (88.1), (94), (97), and (113), there follows the complete expression for /V 3 



a 2 u-^ 2 N, =| [(1 + C, - C2)- 1 / 2 +(1 - Ci- C 2 )- ll2 ](l - S)" 1 ' 2 



77/2 



Vl+2£ 



+ Xo 



+ - [(1 + C, - C 2 )" 1/2 — (1 — C — C 2 )- 1/2 ](l - S)- 1 ' 2 



77/2 



Vl-h2£ 



+ Xi 



(114.1) 
(114.2) 



3 3 

+ (<// + 77/2)^3 —^dC 2 Q cos 1// + — C|<? 2 sin 2i// + periodic terms of 0(/i). (115) 
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Here B,i is given by (95) and x<> and Xi are given by (114). 

To check this expression for 8 = 0, note that then Ci = 0, £==0, S = r)o, u = (a 2 > — a 2 ) -1 a 2 ?? 2 ), 
so that a 2 u~ 112 becomes (a 2 , — a 2 ) 1 / 2 ^ 1 . We have to compare with eq (6.64) of V1961b. Now we 
have already checked #3 and the periodic terms in (115). The check will thus consist only in 
showing that the terms involving xo and Xi in (1 15) reduce to the expression (1 — 7?o)~ 1/2 (1 — T7^ 2 ) _1/2 X 
of V1961b. 

Since Ci = for 8 = 0, the term involving xi disappears, because of the minus sign inside the 
brackets. Since C 2 reduces to 77^ 2 , the first line becomes 



(i-T,r i/2 (i-T) 2 - 2 )- i/2 (|+X")- 



Now the 7r/2 is an additive constant, arising only because we now have the lower limit 171 in the 
^-integrals instead of zero. We shall later get rid of it on replacing the constant /3' 2 by another 
constant /3 2 . Thus we now have only to verify that Xo reduces to x when 8 = 0. When 8 = 0, 
e 2 = (l — P) l Q reduces to 7?o, and we find from (114.1) and (104), with £ = 0, that 



cos- tp— (1-17D) sin- iff 

cos 2xo = t , ■ ., . sin 2xo 

l-m sin- \\f 



2VI — r] 2 ) sin \\t cos iff 
1 — 175 sin 2 1// 



But the x of V1961b, by p. 516 of V1961c, satisfies 



cos x : 



cos l// 



VI — 7)1 sin 2 1// 



smx= 



Vl — T)~ sin \\f 
V 1 — rjt) sin 2 1// 



One verifies easily that cos 2x<> = cos 2x and sin 2x<)=sin 2x« Thus in the limit 8 = 0, the quantity 
Xo and thus the integral Nw both reduce correctly. 



7. Assembly of Equations for the Inversion 

From (11), (67), (78), and (115) of the present paper and (5.30), (5.35), and (5.60) of V 1961b, we 
have 

t + fii= (— 2ai)" l/2 [biE + a(E — e sin E)+A x v + A\\ sini;-h^i2 sin 2v] 



+ C 2 Q!.7 1 tt 1 / 2 



1// + -) #,' + £,, cosi|r + fii2 sin2i// + #,, cos3i|r + Bi 4 sin 4i// 



(116.1) 



f3!>a 2 ] =-(-2a l )- l l 2 



^here 



A>iV-\- V A'zk sin kv 



+ ct7 1 a 1 / 2 



TT 



1// + -) B2+B21 cosi/> 

+ 5 2 2 sin 2^-rB 23 cos3^ + B 2 4 sin4i// 



16.2) 



Bu = 2PQ-^C x Q\ 
Bi 3 =~C 1 ^ 3 /24, 



*..— t$+\w 



fii4 = C2^/64 
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(116.3) 



S 2 3 = -j|c 1 C 2 (? 3 , B 24 = ^Ci0 4 (116.4) 

and where B' 2 and B 2 are given by (76) and (65). 

Note that by V1962 the restrictions on p. 176 of V1961b, concerning orbital inclinations, have 
been removed. 

Before proceeding further, it is convenient to eliminate the terms above involving nj2. If we 
lump the term — c 2 *^ 1 ^ 2 Z?,'7r/2 into )8i, the physical meaning of /3i will not be essentially changed. 
In the limit ,/ 2 — > 0, this change will vanish, since c 2 = 0(J 2 ), and — fit will still reduce in the elliptic 
limit to r, the time of passage through perigee. Also, for J 2 = and /3 = 0, u— 1 and B 2 = 1, so 
that the term — a 2 x u x ' 2 B 2 rrl2 reduces to — 2 tt(x 2 x . Since /3 2 — » w + tt/2 in the elliptic limit, if we 
put 

P 2 = P' 2 -u xIi B 2 ttI2, (117) 

the new constant (3 2 will reduce in the elliptic limit to co, the argument of perigee. 
With these changes, we then have 

t + /3, =(- 2a x )- x l 2 \b x E+ a(E- e sin E) +A&+ A u sin v + A ]2 sin 2v\ 

+ c 2 a L 7 1 ^ 1 / 2 [^;i// + ^,,cosi// + /? 1 .>sin2i// + 5 l:5 c<)s3i//-f B, 4 sin4<M (118.1) 

4 
/^c^ 1 =-(- 2a)- 1 / 2 [^ 2 ^ + 2 ^ 2/ , sin At;] 

f a^i^&^ + fti cos 1// + B22 sin 2i// + fi 23 cos 3i// + B 24 sin 4t//]. (118.2) 
To carrv out the inversion, we write 

E = M S + E l) , v = M s + E Ih i/> = t// s + ^„ (119) 

where the subscript 5 means "secular" and the subscript /> "periodic." 

8. The Secular Solution 

To obtain the secular terms M H and i//. s , we put E = v= z M s in (118) and discard all the sines 
and cosines, thereby obtaining for M s and \fj s the pair of linear algebraic equations: 

{-2a x ) x lHa + b x ^ AMU^ c 2 a- x a^B\^ = t + ^ (120.1) 

-(-2ai) 1 l 2 A 2 M g + a 2 1 l 2 u 1 > 2 B 2 ijJs = /3 2 a 2 1 . (120.2) 



Their solution is 



^here 



M 8 = l +2irvit (121.1), i// s = /o-f^, + 27r^, (121.2) 

27r^ 1 = (-2a 1 ) 1 / 2 (r/ + /; 1 +^ 1 -hc 2 ^ 2J e;^ 2 - 1 )- 1 ( i 2 2.1) 
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27ri/ 2 = a 2 M - ll2 A 2 B* l (a + /;,+/!, + c 2 A 2 B[Bz 1 )- 1 
k = 27Tv^-c^ 2 a7^B\B7^) 
lo+go = 2irv 2 [pi+P2at 1 (a + bi + A l )Az 1 ]. 

9. The Fundamental Frequencies 

By V1961a and pp. 188 and 189 of V1961b, the above v 9 s should satisfy 






(124.1) 



()(X\ 
' $2 



(122.2) 
(123.1) 
(123.2) 



(124.2) 



where j\ and j 2 are the action variables corresponding to p and 77, respectively. From VI 961 h, 
they should thus satisfy 



Vl= OW22-./V21) l J*22, 



(jl\J22-jl2J2l) l jl 



(125) 



y n = 27r(-2a: 1 )- 1/2 (a + ^.+^.), y, 2 = -27r« 2 (-2a,) ""M, (126) 

;ti = 2cWi(i}o), 7*22 = 20:2^2(170). (127) 

Now when 77 = 170, we have i// = 7r/2, by (45.1), (47), and (48). Then, by (71), (78), and (127) 

j 2] = 27rc 2 a 2 ] u^ 2 B f u j 22 = 2nu ll2 B 2 . (128) 

By means of (125) through (128), we readily verify that the fundamental frequencies c)a\/f)j\ and 
c)a 2 /dj 2 agree with the expressions given in (122). 

10. Periodic Terms 

In calculating the periodic terms through the second order, we first set 

E p = £0 + £1 + £2, v p = vo + wi + 1> 2 , ^ = i/*o + </*i + ^2, ( 129) 

where the subscript denotes the order of a term in J 1. We shall calculate the zeroth, first, and 
second order terms in succession. 

Before we can do so, however, we need a table of the order of the various coefficients and con- 
stant quantities. Including the quantities necessary for the right ascension ip also, we find 



Quantity 


Order 


Quantity 


Order 


A 2 


a 


H, 


./■> 


B\ 


A 


in 


J\ 


H, 


A 


A u 


A 


0,2 


A 


ti. 


A 


6, 


Jl 


A 24 


A 


/. 


Jz 


#21 


A 


/>,. 


Jl 


B' 23 


A 


A 2i 


h 


B 24 


A 


■\,,_ 


h 






fil3 


h 






Hu 


h 






II,, 


Jl 
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10.1. Zeroth Order Periodic Terms 

To find £o, v () , and i// () , place E = M S + E {) , v = M s + Vo, \p = \jj 8 + ifo in (118), retain secular terms 
with full accuracy, and discard all periodic quantities of order higher than 75. However, as 
cautioned in V1961b, keep all terms that are retained to as high an accuracy as possible; such a 
procedure simplifies the resulting formulas. Then 

f + fr^^aO-^Ka+fei)^ (129.1) 

j8 2 a2 1 =-(2ai)- 1/ M 2 (M s + t;o)4-a:2 1 " 1/2 fi2(^ + ^o). (129.2) 

Subtraction of (120.1) from (129.1) then gives 

(a + bi)E - ae sin (M s + E ) = 0. (130) 

On dividing (120) by a + &i and adding M s to both sides of the result, we find a Kepler equation for 
M s + Eo, viz, 

Ms + Eo- e' sin (Ms + E<>) = M s (131.1) 



ae 



<e<l. (131.2) 



a~\- b\ 
Solve this Kepler equation for E = M s + Eq. Then find v = M s + Vq by the anomaly connections 

cosE — e . (1 — e 2 ) 1/2 sin E 

cos v = - = sin v = — =— • (132) 

1 — e cos L 1 — e cos E 

To find i// , subtract (120.2) from (129.2) and solve for t// () . The result is 

i//o = a 2 (- 2a x ) -^u-^A 2 B- l v». (133) 

As in V1961b, this gives i// = vo -h OU2), but it is important to retain the coefficient with full 
accuracy, for later simplifications. 



10.2. First Order Periodic Terms 

We now place E = M s + E + £1 , 1; = Af, + 1; + f>i , 1// = i// v + t// + 1//1 in (118), discarding periodic 
terms of order 7|, and use our knowledge of M s , i// s , E , v , and i// () to find £1, i>i, and ipu From 
(118) we obtain 

t + p x =(-2a x y^ [(a-r-ix) (M s + £ + £1) - ae sin (flf f + £ + £i) 

+ ^(M, + ^ )]+c 2 a L 7V/ 2 [i?J(^H-^o)+Bi2 sin (2i//, + 2t/> )] (134.1) 

ftas"^- (-2a,)- 1 / 2 [^(Af.+ vo + fi) + ^21 sin (Af« + i;„) + /4 22 sin (2M, + 2t/ )] 

+ a 2 " 1 a 1 / 2 [&(*. + 0b + M+&2 sin (2i//, + 2i// () )]. (134.2) 
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On subtracting (120.1) from (134.1), transposing some terms, and dividing through by (— 2«i) 1/2 
(a + /;i), we find: 

E,^ E } - e' sin(M s + Eo + £,) = - (a^ b.r'lAm^ c 2 a 2 H-2a l ) 1 l 2 u^ 2 B[ijj {) 

+ c 2 oq\- 2a,yl 2 a l l 2 B 12 sin (2i// 8 -f 2i// () )], (135) 

where e' is given by (131.2). By (133) the terms in vo and i// consolidate into (A x + c 2 A 2 B[B 2 x )v$ 
and (135) becomes 

M s + £ + £'i-e , sin(M s + £o + £i) = M s -hM 1 , (136.1) 



M 1 =-(^ + />,)-HU.+c 2 ^25;5^M^)+c 2 a7^-2a 1 ) 1 /^/ 1 /^ 12 sin(2i//. s . + 2i/;o)J. (136.2) 

By p. 194 of V1961b, the solution for £i is then given by 

r Mi e' /tffsin(M,-hE () ) 



" J 1 -e' cos (A/,. +£ ) 2 [1-e' cos (M 8 + E )f* W?) ' (137) 

To find i;i, insert the known values of E = M s -\-Eo~\- E \ and e into the anomaly connections (132) 
and solve for v = M s + l>o + V\. 

To find i//i, subtract (120.2) from (134.2). The result is 

a2 1 H 1/2 [fl 2 W*> + M + B22 s^ (2i/; s .+ 2i// () )| = (-2a 1 )- 1 / 2 [^ 2 (z; ( )-hi;,) 

+ ^2i sin (M 8 + vo)+A 22 sin (2M , + 2v )]. (138) 
By virtue of (133) the terms proportional to i// () and Vo cancel out. On solving for i//i, we then find 
iPx=-B 22 B- 1 sin (2i//., + 2i// () ) + a 2 (-2a 1 )- 1 / 2 w- 1 / 2 ^- 1 [^2^i 

+A 2 i sin {M 8 + v )+A 2 2 sin (2M* + 2i>o)]. (139) 

10.3. Second Order Periodic Terms 

To find £ 2 , t> 2 , and i// 2 , insert E=M 8 + £ + £1 + £2, v = % + t>o + i>i + v 2 , and 1//= t/f* + *//<> + 1//1 + i/i 2 
into (118), retaining all periodic terms of order J| and discarding all those of higher order. We 
obtain 

t + p 1= (-2a l )- l ' 2 [(a + bi) (M s + E {) + Ei + £2)- ae sin (M 8 + £ + £1 + £2) 

+ ^i(M« + v +t;i) + ^ii sin (M g + i;o) + ^i2 sin (2M s -h2i; )] 
H-c 2 a 2 - 1 M 1/2 [B;(^ + ^o + W + Bu cos (i// s + «|r ) + Si2 sin (2^+2^ + 2^i) 

+ B13 cos (3^,+ 3*b) + Bi4 sin (4i//, + 4i// )]. (140.1) 
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fca^ =- (-2a ] )- l t 2 [A 2 (M s + v {) + vi + v 2 ) + A 2l sin (M s + v + Vi) + A 2 .> sin (2M S + 2v + 2t;i) 

+ ^ 23 sin (3M s 4-3t;o)-f^24 sin (4Af„ + 4t;«)] + a 2 - 1 a 1/2 [fi 2 (^ + i/fo + 4>i + <fe) + S 2 i cos (ife+i/fo) 

+ B 22 sin (2i//, + 2i// + 2i// 1 ) + J e 2 : J cos (3ifc + 3i// ) + B 24 sin (4<|/, + 4i// () ) |. (140.2) 

To solve for E>, subtract (134.1) from (140.1). We find 

(-2«i)- ,/2 [(a + b x )E 2 - ae sin (M s + E + E x + E 2 ) + ae sin (Af s -f-£ + £i)+^it>i+^ii sin (Af. + Vo) 

+A X2 sin (2M s .-h2i;o)j + c 2 a 2 ~ 1 a 1/2 [J5;i//,4- J e il cos (^-f^o) + Bi 2 sin (2i//, + 2i// + 2<//,) 

-B 12 sin (2i|f s + 2iM + #i3 cos (3i//, -h 3i// ) + ^i4 sin (4i//,. + 4i//<>)]. (141) 
Now 

sin (M« + E« + £,+£*) -sin (Af, + £ + £i) = £2 cos {M s + Eq + E x ) + 0(J$) (142.1) 

and 

c 2 sin (2i// s H-2^o + 2i^,)-c 2 sin (2i//, + 2i// ) = 2c 2 i//, cos (2i/> s + 2iM + WI), (142.2) 

where £ 2 and i//i are in radians. On inserting (142) into (141) and solving for E 2 , we find 

* 2= W cos (*#. + £„ + £,) ( 143J ) 

M 2 =-(a+6,)- 1 [^,?;,+^n sin (M,. + z; ) + ^, 2 sin (2Af, + 2t*>) 

+ c 2 a L 7M-2a,) 1/2 // 1/2 {#,>,+#,, cos (^ + W + 2Bi2^i cos ( 2 1//., -h 2 1// ) 

+ #,, cos (3i//, + 3i//„) + #i4 sin (4i//, + 4i// () ) } ]. (143.2) 

Comparison with (8.45) of V 1961b shows that a term in sin ( 2 i//. s - -h 2 1//< > ) is now missing. The 
reason is that, by using the full accuracy for B\ 2 in the first order periodic terms, we have already 
included it in M\ . 

To find v 2 , insert E = M s H- £o + E\ + E 2 in the anomaly connections (132) and solve for v 

= Af, + t;o + t;i + !;2. 

To find i|/ 2 , subtract (134.2) from (140.2) to obtain 

-(-2ai)-^ 2 [A 2 v 2 -\-A 2 i sm{M s +v {) + Vi)-A 2 i sin (M, + t*>) 

+ ^a sin (IMs + 2«o + 2vi) - ^22 sin (2 Mi, + 2v ) 

+ /i 23 sin (3M, + 3v ) + ^24 sin (4M,. + 4v ) I 

+ a.7 1 // ,/2 [^i//2 + ^2i cos (i//.v + i//o) + #22 sin (2i//., + 2i//«, + 2i//,) 

- #22 sin (2i// v + 2i// () ) + B 2:i cos (3i//. v -f 3t//<>) + B 2A cos (4i//. s . + 4i|/ () ) ] = 0. (144) 
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On inserting 



/•j, sin (M 8 + t*, + vi) - A n sin (M, + 1; ) = A, x ih cos (/^, + v n ) + 0(JD 



A22 sin (2/1/., + 2v + 2ri ) - ^22 sin (2M. S - + 2v ) = 2A 22 v ] cos (2M, + 2t*>) + 0(7| ) 
B n sin (2i//. s .+ 2i// + 2i/i 1 )-/i 22 sin (2i^.,-H 2i// ) = 2fi 22 ^i cos (2t// g + 2i/f ) + 0(7 2 3 ) 
into (144) and solving for 1//2, we find 

i// 2 = - Ba'fBai cos (i//.s. + i//,)) + 2Z? 22 i//, cos (2i//,. + 2i//„) + B 23 cos (3i// x + 3i//„) + #24 sin (4i//, + 4i// ) | 
+ a 2 w- 1/2 (-2a 1 )- 1 /2^ 2 - 1 [A 2 v 2 + A 2 iVi cos (M,. + *; () ) + 2/t^i cos (2M, + 2i; ) 

+ A 2:i sin (3M S + 3i> ) + /* 24 sin (4M* + 4z; () ) |. (145) 
This completes the solution for the E, v, and 1// and thus for the spheroidal coordinates p and 77. 

1 1 . The Right Ascension 

To find if we have to insert (115) for /V 3 and (5.60) of V 1961b for R A into (11.3). Then 

r 4 

<p = Ps-<?cta(—2ai) ,/2 Azv + ^Asj sin jv 

L j =1 J 

+ a 3 a2 1 w 1/2 j(l-S)- 1/2 (A, + A 2 ){|(l + 2£) l/2 + X ] 

+ (A,-/i 2 )|-|(l + 2£) 1/2 + x } +(|+i//)B 3 -|c,C 2 <?cosi//-f^q(? 2 sin2^j- (146) 

Here /*;} is given by (5.61) through (5.65) in \ 1961b. Because of typographical errors there, 1 
include here the correct expressions for /:>_> and ^33: 



4w = (l-e 2 ) 1/, P :i 



e 1 3 6, 2 _ 2 A? 4 , 3 2 \/6? 



e 3 ^ e 3 (b% , 9 



12/; 3//\2 



(146.1) 



Here b\ is simply the B in (34.1) of this paper. 
The symbols h\ and h 2 are abbreviations: 

1 



/ ?1 -^(l-fC 1 -C2)- ,/2 



147.1) A 2 ^^a-Ci-c 2 )- i/2 - 



(147.2) 



By (18), (21.2), (26). and (31), we find for a, 



a3 = [Atpo(l -S)| 1/2 sgna.-j 



1- 



c 2 S \ poj \ aopo/ 



(lo/h) 



aopo 



(1-25) 



1/2. 



The error here is of order J I. With the same error, we can simplify this to 



a 3 = [/i,po(l — S)] 1/2 sgna 3 



aopo \/W \ (ip ap 



1 2 



(148) 



(149) 
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Alternatively, when ai is known, we may write 



S\i/2 
a 3 = a 2 \ 1 — ^1 sgn a 3 ■ 



(149.1) 



In (146) the constant /3 3 would reduce to (1 — 2 77 sgn oi 3 in the Keplerian case, as one may see 
by solving the two body problem with lower limit T7i(= — 170) in the 77-integrals. Also in (146), 
consider the constant term: 

I ttw^ l u^(\ - S)- 1 ' 2 (1 + 20" 1/2 [/n + h 2 ~ (/m - k*)] + I 77# 3 . 

For J> = and ,/ 3 =0, we have w = l, £ = 0, h\ J th 1 = 1, hi — h 2 = 0, Z? 3 = 0, 5 = sin 2 /, and asaj 1 
= cos /. Thus it then becomes 2 77 cos /| sec /| = 2 77 sgn a 3 . If we lump it into the constant 
/3 3 and call the sum /3 3 , we then find that /3 3 reduces, in fhe Keplerian case, to ft, the right ascension 
of the node. 



Our amended equation thus takes the form 



<p = fe-c 2 a 3 (-2a 1 )- 1 l 2 



A 3 v+ ^Azj s'mjv +a 3 a 2 - 1 ^ 2 Ul-S)- 1 l 2 [(h 1 + h 2 )x o 



+ (hi-h 2 )x 1 ] + B3ilf~ C x CzQ cos 1// -h ^ C|<? 2 sin 2i//j, (150) 

where x o an d Xj are given by (114) and E' 2 and E' s in (114) by (104), with e 2 = Q(l — P) x 
and6 3 = <?(H-P)- 1 . 

It is not difficult to verify that (150) reduces to (8.50) of V1961b in the limit 8 = 0. A much 
more delicate check is to show that (150) makes <p a constant in a polar orbit, with jumps of 77 as the 
satellite goes over a pole. 

As an orbit approaches being a polar orbit, a 3 — »0 and S — > 1, so that 



<p -* fc + otaoq l u l l*[(ht + h 2 )x + (k ~ Wxjd- S)~ 1/2 . 
Now by (148), since (/xp ) 1/2 = o:2, we have for S— > 1 



(151) 



a^a l u ll2 (l — S) 1/2 =sgn a 3 



c 2 (25/po) 2 



dopo c l 



dopo 



1/2 



(152) 



since by (31) a — > 1 as S — » 1. Also, for S= 1, we have by (32) that (152) becomes 

asa^a^U-S)- 1 ^ sgn a 3 [(l-C2) 2 -Cf] 1/2 a-C 2 )- 1/2 . (153) 

To simplify the term in (151) involving the x's, note that by (102.2), (32.1), and (32.2) 

2£ = (l-C,r 1 G. (154) 

Then, by (147) and (154) 

(1 - 2£)~ 1/2 = 2W1 - G) 1/2 (155.1), (1 + 2Q^ = 2*i(l - C^i\ (155.2) 
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In (114) let us now abbreviate El,(\\i + it 12) to E! 2 and E!^- tt/2) to £3. Equations (114) and (155) 
then yield 



X=a-C>)n^>E^h x E[,) 
X =a-C 2 ) ll2 (h i E! 2 -h l E' 3 ), 



from which it follows that 

(*i + ^)X + (*i - Wx, = 2(1 - C 2 ) ll2 h l h2(E! 2 + £ 3 ) 



(156.1) 
(156.2) 

(157) 
Hl-Q^Kl-C^-qi-^f' + ^a), (158) 



after insertion into (157) of the expressions (147) for h\ and /i 2 . On inserting (153) and (158) into 
(151), we then find that for a polar orbit 



(p-*fc+ 5(E! 2 + E' z )sgn a- s . 



(159) 



To discuss this expression, note that e2 = Q(\— P) ' and en = Q(l + P)~ l are both unity for a 
polar orbit, since by (32.3) P = when S=] and since by (47) Q then equals S=l. 

To see how E' 2 -\- E' :l behaves as S —> 1, we have to discuss the function E'(f), where the argu- 
ment /plays the role of true anomaly and the function E' the role of eccentric anomaly. 

Now by their definitions E' and/agree in value at all multiples of 77. Also 



dE' Vl-e 2 



df \ + e cos/ 



Tin 



dE' (\ 



1/2 



df \\ + e) 

1 + e V 
1-e 



>0ase->l, /=2t7t 



00 ase-^1, /=(2T-r-l)7r 



The plot off versus /for e=l is thus as follows 



(160) 



(161.1) 



(161.2) 



FIGURE 1. Plot of eccentric anomaly versus true anomaly 
for e =/. 



7 / 



77 2TT 37T 47T 5 77" 67T 77T 



43 



For e=l, one sees that E'{f) remains constant from /= (27 + 1)77 to (27 + 3)77, but then jumps 
by 2?7. 

77 77 

Let us assume sgn a^ = \ and consider what happens as t// increases from — — + e' to — — e'. 

In this range *//+ tt/2 increases from e' to 77 — e' and i// — 7i/2 from —77 + 6' to — e'. In this range, 
therefore, both £"(1// + 77/2) and E'{ijj — tt/2) remain constant, so that <£> remains constant, by (159). 
As 1// increases beyond 7r/2, however, i//-h 77/2 passes through an odd multiple of 77 and iff — tt/2 
through an even multiple of 77. Thus £"(i// + tt/2) jumps by 277 and E'(\fj — 77/2) remains constant, 
so that, by (159), <p jumps by 77. This is at the north pole. A similar jump takes place at the south 
pole. These results are of course trivial for a polar orbit, but they constitute an important test 
that the equation (150) for the right ascension had to pass. 



12. The Final Algorithm 

To summarize results, assume that /*,, r e , J2, J3, and the orbital elements a, e, S, /o, go, and ^83 
are given. (Here 5 is an element that would reduce to sin 2 / — T7o of V1961b, but with J 3 incorpo- 
rated into the reference orbit it may go negative for orbits close to equatorial.) 

To calculate the reference orbit, compute 

c 2 = rjj 2 (l-|ji./A 8=-| rJ^J 3 > for J 3 < p = a(l-e 2 ), 

A by (40.1), and B by (40.2). Then 

«o = a - §A, p = <h l (B + ap - 2Aa - c 2 ), a 2 = (/W /2 , u by (31), C 2 by (32.2), C, by (31.1), P by (32.3), 
b\ = — 2 A, b 2 = B 112 , ol\ =— I M-a ( 7'. Then compute 

A u A 2 , As, A n , A 12 , A iu A 22 , A 23 , A 24 , A 3U A 32 , A 33 , A 34 by p. 198 of V1961b. 

Now compute Q = (i^ + S) 1 ' 2 , B> by (65), B[ by (76), B, by (95), £ by (102.2), 

e 2 = Q(l-P)- 1 e^QH + P)' 1 
2irv x by (122.1), 2rrv 2 by (122.2), and Z and go by (123). Then 

M s = A, + 2771^, \jj s = /o + g {) + 2irv 2 t 

e '-- 2£-<«<l. 

a + /)] 



Find M s + £"() by solving 



Ms + E -e' sin (M, + E ) = M s . 



Solve (132) for t; = M s + i; and (133) for i// . Then find M, by (136.2) and E x by (137) and solve 
(132)fort; = M,+ t*> + t?i. Find i/>i by (139) and then M 2 by (143.2) and E 2 by (143.1). 

Solve (132) for v = M s + v + Vi + v 2 and find \Jj 2 by (145). Then 

E=Ms+Eo + Ei+E 2 ,v = M s + vo + vi + v 2 , ^=^ + ^ + ^1 + ^2. 
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The spheroidal coordinates p and 77 are then given hy 

p = a(l — e cos E) = (1H- e cos i;) _1 p 

7) = P + Q sin 1//. 

To find the right ascension <p, note first the correct formulas for ^32 and ^33 in (146.1) of this 
paper or on p. 198 of V1961b. Then compute hi and h 2 by (147) and a.? hy (148). Define Zi^y) by 

„,. x e L > + cos y „„ . (1 — e|) ,/2 sin y 

cos E 2 (y) = — — sin £/ 2 (y) = 

1 H- e2 cos y 1 + e 2 cos y 

and the requirement that £3' and y shall be equal whenever y is a multiple of 77. Define E' z (y) in the 
same way, with e% replacing cn. 

Then calculate \ and x, by (114) and if by (150). Having p. 77, and ^>, we may now find the 
rectangular and spheroidal coordinates by means of (1.1) and (1.2). 

13. Further Remarks 

The present inclusion of J3 into the separable orbit has apparently led to new secular terms 

13 3 

of order J\. These terms, — - C\P-\- — C\Q- in 62 and — — C\ in 63, are proportional to J%l.r 

Z 10 o 

= O(J'i) and lead to terms in l/. s , i//. s , and <p. Since no previous treatment of third harmonic perturba- 
tions has led to such secular terms, the result requires checking. 
To do so, note that by (2) and (3) the present c 2 is given by 

c 2 = cf-8 2 

where Co = r 2 J 2 . It turns out that the 8 2 terms implicit in the new c 2 are just sufficient to cancel 
the new terms in § 2 . 

Since c 4 = c,1 — 2c 2 S 2 -\r 8 4 , terms of order c 4 cannot cancel terms in 8 2 . In expanding the fun- 
damental frequencies it is thus sufficient to carry terms only through c 2 and 8 2 , in checking this 
point. A rather long calculation then shows that 

M s = 2Trvi = n + O(c 4 ) 



lfj s =27TV2— Hi) 



(p = 2ttv;i = no sgn a. 



1+ !4< 4 - 5S » 

4 p z 

3 c 2 tM _ _ 3 c; 



+ 0(c 4 ) 



l+-^(4-5S)-2^d-S)^ 



4- 0(c 4 ). 



Here n {) = p, ll2 a^ 312 , the mean motion corresponding to ao = — fx{2a\) \ where ot\ is the total energy. 
This means that, if we solve the problem for given initial conditions without using the third harmonic 
and then re-solve it with inclusion of the third harmonic, there will arise no secular terms of order 
J 2 j7 2 , provided that we use the same initial total energy. In turn, this means that if we use the 
same initial coordinates, we must change the initial velocities slightly or that if we use the same 
initial velocities we must change the initial coordinates slightly. 

It does not appear worth while to do the extra calculation that would give us final formulas 
manifestly free of secular terms of order JiJ* 1 - If we use the present results, we shall obtain secu- 
lar terms correct through order J\\ this limitation comes from the error of order J\ in (31.2). To 
obtain this secular accuracy, we find from (122) and (150) that we should keep terms through order 
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J\ in B 2 and B 3 and through order )\ in B[. By improving the solution of the cubic (27), we could 
calculate secular terms to any order of accuracy. Periodic terms are correct through order J 2 ; 
they also may be carried to a higher order if the need ever arises. 

To obtain the orbital elements one may proceed as follows. From initial values determine 
«i, a 2 , and a 3 by the appropriate modifications of eqs (3.8), (3.9), and (3.10) of V1961. They are 

«i = 2 u\ — n(pi + rji8) (p 2 + c 2 ^ 2 )" 1 

a 3 = r 2 cos 2 Oi(pi = %\ji — yixi 

a 2 = (1 - 7)\Y l M + c 2 i) 2 ) 2 i) 2 + a 2 - (1 - n 2 ) (2aic 2 7jf + 2/xi|f6)]. 

Here the equation for a 2 follows from (10.2), (13.2), (50), and (59.1) of V1959. The subscript i means 
initial value and u\ is the initial speed. 

From the a's one can determine a , eo, po, and to from eqs (21) and (28) of the present paper 
and then determine a, e, and 170 = sin / by means of the equations on pp. (172) and (173) of V1961b. 
Knowing a, e, and / to this approximation, one can then determine the values of j8i, ^2, and /3 3 or 
of /o, go, and /3 3 by fitting the values of p, 17, and <p at various times, still using the equations in 
V1961b. 

To improve these elements, now calculate the orbit, with elements a, e, S = sin 2 /, /o, go, and 
^83 by means of the present paper. Observe the residuals and use them, with the aid of least 
squares, to obtain a differential correction to improve the orbital elements. Here a word of cau- 
tion is necessary. If the orbit is close to equatorial, S may go negative. This is all right, since 
it is only Q 2 = P 2 + S that cannot become negative, with the present method. 

Then repeat the process, until the elements converge to fixed values. One or two iterations 
should suffice, according to Bonavito (1964). 

The present method of accounting for J 3 ought to be an improvement over the perturbation 
theory of V1963, for two reasons. There appear no small denominators, such as those containing 
e or sin / in the perturbation theory. All periodic terms are correct through order J\. 
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